The amorphous sohd state: a locally stable thermodynamic phase of randomly 

constrained systems 



Horacio E. Castillo^"' Paul M. Goldbart^'') and Annette Zippelius^^^ 
(") CNRS-Lahoratoire de Physique Theorique de I'Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris CEDEX 05, 

France 

Department of Physics, University of Illinois at Urbana- Champaign, 1110 West Green Street, Urbana, IL 61801-3080, USA 
Institut fiir Theoretische Physik, Georg August Universitdt, D 37073 Gottingen, Germany 

(May 21, 1999) 

The question of the local stabihty of the (replica-symmetric) amorphous solid state is addressed 
for a class of systems undergoing a continuous liquid to amorphous- solid phase transition driven 
by the effect of random constraints. The Hessian matrix, associated with infinitesimal fluctuations 
around the stationary point corresponding to the amorphous solid state, is obtained. The eigenvalues 
of this Hessian matrix are all shown to be strictly positive near the transition, except for one — 
the zero mode associated with the spontaneously broken continuous translational symmetry of the 
system. Thus the local stability of the amorphous solid state is established. 



I. INTRODUCTION 



In recent years, a theoretical approach has been de- 
veloped for the problem of the liquid-amorphous solid 
phase transition in systems of randomly crosslinked flex- 
ible linear macromolecules [0-^. This approach starts 
from a semi-microscopic model for the macromolecules, 
and takes into account explicitly both the thermal fluctu- 
ations at nonzero temperature and the quenched disorder 
due to the random nature of the crosslinking. It is based 
on the Deam-Edwards formulation of the statistical me- 
chanics of polymer networks and borrows some con- 
cepts and techniques that have in the past been used to 
study the problem of spin glasses Q . In the framework of 
this approach, it is possible to define an order parameter 
that probes random static particle-density fluctuations, 
and is thus able to detect the transition between the liq- 
uid and the amorphous sohd state 

Along the way, it has been recognized that there ex- 
ists a class of systems (including, e.g., end- linked flexible 
and stiff polymers , and crosslinked higher-dimensional 
manifolds |^) that display identical transitions, and a 
general Landau theory was formulated to describe this 
transition Q. This Landau theory was constructed us- 
ing only symmetry considerations and the assumption 
that the phase transition should be continuous (and thus 
dominated by the long-distance behavior of the system). 

By expanding any of the semi-microscopically de- 
rived free energy functionals corresponding to the above- 
mentioned physical systems in powers of the order pa- 
rameter and gradients, one recovers this general Landau 
free energy functional, up to terms that play no role in the 
mean field theory, in the vicinity of the transition. In the 
study of fluctuations, as we shall show later, these addi- 
tional terms do not alter the physical picture. The slight 
difference between the general Landau theory and the 
microscopically-derived theories is due to the fact that. 



whereas the former theory only allows states that have a 
disorder-averaged particle-density that is spatially homo- 
geneous, in the latter states with spatial inhomogeneities 
of the density are in principle allowed, but are ultimately 
suppressed by the repulsive interparticle interactions. 

Let us begin by summarizing the results of the mean 
field theory of the liquid-amorphous-solid transition 
driven by random constraints: (i) for densities of con- 
straints smaller than a critical value the system is in the 
liquid state and all particles are delocalized; (ii) for densi- 
ties of constraints larger than the critical value the system 
is in an amorphous solid state, characterized by emerg- 
ing random static density-fluctuations; (iii) at the critical 
density of constraints there is a continuous phase transi- 
tion between the liquid and the amorphous solid states; 
(iv) in the amorphous solid state a positive fraction of 
the particles is localized around random mean positions 
and with random r.m.s. displacements; (v) in the amor- 
phous solid state, close to the transition, the fraction 
of localized particles is proportional to the excess of the 
crosslink density beyond its critical value, and the typi- 
cal localization length diverges at the transition like the 
excess crosslink density to the power —1/2; (vi) when 
scaled by the mean value, the statistical distribution of 
localization lengths is universal near the transition and, 
consequently, the dependence of the order parameter on 
the wavevectors also has a universal scaling form. 

Although in the amorphous solid state translational in- 
variance and rotational invariance are broken at the mi- 
croscopic level, because a fraction of the monomers are 
localized, the average density is uniform, and the system 
is macroscopically translationally (and rotationally) in- 
variant (MTI) The liquid state is, of course, not only 
macroscopically but also microscopically translationally 
and rotationally invariant, as each individual monomer 
density is uniform over the container. Moreover, as in 
the case of some of the mean field solutions encountered 
in spin glass systems, both the liquid and the amorphous 
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solid states are replica-symmetric |Q. However, as will 
be shown later, the symmetry of these states is even 
larger, because the replica symmetry combines with the 
rotational symmetry (in d-dimensional space) to produce 
rotational symmetry in a nd-dimensional (i.e., replicated) 
space. 

It is well known that in the spin glass case the 
early replica-symmetric spin glass solution fictl for the 
Sherrington-Kirkpatrick model was later found to be lo- 
cally unstable at temperatures below the critical temper- 
ature by de Almeida and Thouless and was super- 
seded by the replica-symmetry-breaking solution jl^ dis- 
covered by Parisi, with its elegant interpretation in terms 
of symmetry- unrelated pure equilibrium states . This 
replica-symmetry-breaking solution has a nonnegative 
entr opy at all temperatures and is locally marginally sta- 
bleQ. 

By analogy with the spin glass case, and consider- 
ing also the fact that random topological constraints 
can, in principle, lead to a partitioning of the configura- 
tion space of three-dimensional macromolecular systems 
into ergodic regions that are not connected by symme- 
try operations , it would not be entirely surprising if 
in systems undergoing a liquid-amorphous-solid transi- 
tion under the effect of random constraints, the replica- 
symmetric stationary point of the free energy correspond- 
ing to the amorphous solid state turned out to be unsta- 
ble and had to be superseded by a less symmetric solution 
of the mean field equations. 

However, the Deam-Edwards distribution ||,^ used to 
model the disorder favors sets of constraint locations that 
are associated with highly probable configurations of the 
unconstrained system. Thus, although frustration may 
in principle be present, the disorder distribution tends to 
discourage it [ p^ . 

The purpose of this Paper is to show that the mean 
field amorphous solid state is, in fact, locally stable, at 
least near the transition. More specifically, the Hessian 
matrix that describes changes of the free energy func- 
tional for infinitesimal fiuctuations around the stationary 
point corresponding to the amorphous solid state is com- 
puted, and it is shown that, to linear order in the excess 
crosslink density, all its eigenvalues are strictly positive, 
except for a single zero mode associated with the sponta- 
neous breaking of the continuous translational symmetry 
of the system. In order to do this, we construct a descrip- 
tion for the space of fluctuations around the stationary 
point that allows the eigenvalue problem for the Hessian 
matrix in replicated space to be reduced, in essence, to 
an integral eigenvalue equation in one dimension. 

The rest of this Paper is organized as follows: In Sec. || 
we compute the Hessian matrix for the Landau free en- 
ergy functional around the amorphous solid state, and 
make use of the continuous symmetry of the problem to 
identify a basis set that significantly si mplifi es the eigen- 
value equation for the Hessian. In Sec. Ill wc find posi- 
tive lower bounds for all the eigenvalues of the Hessian, 
except for the zero mode that is present due to the spon- 



taneously broken symmetry. In Sec. IV we extend these 



results to allow for fluctuations in the particle density in 
the case of the microscopically-derived free energy func- 
tional for randomly crosslinked macromolecules. Finally, 
in Sec. ^ we present our conclusions. 



II. LANDAU THEORY: HESSIAN MATRIX 

In this section we first summarize briefly the basics of 
the general Landau theory of the liquid-amorphous-solid 
phase transition, as presented in Ref. Q, and later de- 
rive expressions for the Hessian matrix for the free energy 
functional in this theory. 



A. Brief review of the Landau theory 

In a system characterized by static random density 
fluctuations, the appropriate order parameter is lllJ 
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(2.1) 



where N is the total number of particles, (with i — 
1 , . . . , N) is the (d-dimensional) position vector of parti- 
cle i, the wave- vectors k^, k^, • • • , are arbitrary, {■ ■ ■)x 
denotes a thermal average for a particular realization x of 
the disorder, [• • •] represents averaging over the disorder, 
and g is a positive number. 

We make the Deam-Edwards assumption [|j that the 
statistics of the disorder is determined by the correlations 
of the unconstrained system. Under the Deam-Edwards 
assumption, obtaining disorder averages with the replica 
technique amounts to working with the n ^ limit of 
systems of as opposed to n, replicas. The additional 
replica, labeled by a = 0, represents the degrees of free- 
dom of the original system before adding the constraints 
or, equivalently, describes the constraint distribution. 

In the replica formalism, the order parameter takes the 
form 



1 ^ 

i=l 



P 

ri+1 



(2.2) 



Here, hatted vectors denote replicated collections of [d- 
dimensional) vectors, viz., v = (v*^, v^, • • • , v"), their 
scalar product being v ■ w = X]a=o ' (' ' ')n+i 

denotes an average for an effective pure (i.e., disorder- 
free) system of n -I- 1 coupled replicas of the original 
system. We use the terms one-replica sector (Irs) and 
higher-replica sector (hrs) to refer to replicated vectors 
with, respectively, exactly one and more than one replica 
a for which the corresponding vector k" is nonzero. 
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In the Landau theory the order parameter in the 
one-repHca sector represents spatial variations in the 
disorder-averaged mean particle-density, and is always 
assumed to be strictly zero. 

In the stationary-point approximation, the disorder- 
averaged free energy / (per particle and space dimension) 
is given by §,0,0 



with the Landau free energy functional given by 



(2.3) 



2 



^kik2k3 '^l 



(2.4) 



Here e is the control parameter, and is proportional to 
the amount by which the constraint density exceeds its 
value at the transition. The symbol ^ denotes a sum 
over replicated wave vectors k in the higher-replica sec- 
tor. 

For e < 0, the stationary-point equations for the free 
energy functional of Eq. (2.4) only have the solution 



rif, — 5f,Q, corresponding to the liquid state. For e > 0, 
there are two solutions, one (O^ = 5f. g) corresponding 
to the liquid state and a second one corresponding to an 
amorphous solid state, given by 



% = (l-'?)'5L6 + 9^k,o' 
q = 2e/3, 

poo 

uj(k) = / d6l7r(6l)e^'='/2^ 



(2.5) 



Here, the quantity q is the ratio of the number of lo- 
calized particles to the total number of particles (i.e., 
the localized fraction), 7r(^) is a universal scaling func- 
tion that characterizes the distribution p(l/^^) of the 
(inverse square) localization lengths ^ for the localized 
particles, through p(l/^^) — (2/e)7r(2/e^^), and we use 
the definition k = J2a=o ■ "^^^ factor p encodes the 
property of macroscopic translation invariance (MTI) for 
the amorphous solid state, i.e., the fact that the state is 
invariant under common translations of all the replicas, 
or, in more physical terms, that the particles are local- 
ized around randomly located points that have a homoge- 
neous probability of being found anywhere in the volume 
of the system. The scaling function Tr{9) satisfies the 
stationarity condition 

y ^ = (1 - 0) n{0) J' de'n{e'y{9 - 9'), (2.6) 



together with the normalization condition 



d9 7r{0). 



(2.7) 



This normalization condition directly f ollow s from the 
fact that the order parameter of Eq. (2^) has to be 
unity at the origin of replicated wave vector space . 
It is worth noticing that in the limit e — the above 
parametrization of Q reduces continuously to the order 
parameter Qj^ — for the liquid state, as it should. 

Let us now discuss the symmetry properties of the Lan- 
dau free-energy functional. Under independent transla- 
tions of all the replic as, i .e., of — > of + a", the replica 
order parameter, Eq. (|2.2|), transforms as 



k 



(2.8) 



For later reference, let us calculate the change in the or- 
der parameter for the case of small displacements of the 
replicas: 



Sflf^ = fl'j^ — Clf^ = ik-a^f. 



0{a'). 



(2.9) 



Under independent rotations of the replicas, defined by 
cf ^ c'" = i?"cf , and Rij = {R°v°, • • • , i?"v"}, where 
each R" is a rotation matrix in d dimensions, the order 
parameter transforms as 



^k 



(2.10) 



By inserting the transformed order parameter for either 
of t he a bove operations into the free energy functional 
Eq. (2.4), we see that in both cases: 



ndTniin'j^}) ^ ndTn{{nf^}), 



(2.11) 



i.e., that the Landau free energy is invariant under inde- 
pendent translations and rotations of the replicas. 

Anticipating the conclusions of this Paper, that the 
liquid state, which becomes unstable when e is increased 
through zero, is replaced by a stable amorphous solid 
state for e > 0, we now pause to compare some aspects 
of this phase transition with their counterparts in simple 
models of the paramagnct-to-ferromagnet phase transi- 
tions [e.g., the 0{N) symmetric vector (j)'^ model]. We 
shall refer to Fourier components of fields as modes, and 
shall consider the mean-field level of description. At high 
temperatures (for magnetism) and low constraint densi- 
ties (for amorphous solidification) the equilibrium value 
of all modes is zero. As the relevant control parameter 
is changed through its critical value, a band of modes, 
including those of the longest wave length, become lin- 
early unstable, the longer the wave length the stronger 
the instability. In both settings, magnetism and amor- 
phous solidification, stability is recovered by the acqui- 
sition of a nonzero equilibrium value by one or more of 
the modes. For magnetism, there is a zero wave vector 
mode, which is the most unstable mode, and by giving it 
the appropriate nonzero equilibrium value, this mode and 
all others are restabilized (i.e. are no longer wnstable — 
as discussed below, there should and does remain one 
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marginally stable Goldstone mode). For amorphous so- 
lidification, restabilization is more intricate. There is no 
fluctuating zero wave vector mode in the theory to be 
given a nonzero equilibrium value. Instead, the most 
unstable modes have the smallest allowed nonzero wave 
vectors. If these modes become nonzero, as some of them 
do, there is no symmetry-dictated "selection rule" pro- 
hibiting them from acting as "sources" for certain other 
modes, and thus not just one but a large family of modes 
become nonzero, not only including modes that were for- 
merly unstable. (This may be regarded as an analog 
of domain wall formation in wave vector space.) There 
is one further subtlety to the amorphous solidification 
case: to give a nonzero equilibrium value to modes that 
reside in the one replica sector would be extremely en- 
ergetically costly and, from the viewpoint of the Landau 
theory, is ruled out (via an implicit linear constraint on 
the order parameter). This requirement is satisfied by 
giving a nonzero equilibrium value only to modes that 
are MTI, because such modes are prohibited, on symme- 
try grounds, from acting as sources for modes in the one 
replica sector. Thus, stability is restored not by giving 
one unstable mode a nonzero value, and not by giving 
only the unstable modes a nonzero value, but by giv- 
ing a sheet of modes (some unstable and some stable) a 
nonzero value. 



B. Hessian matrix elements 



(2.14) 



Thus we can rewrite the second variation in terms of the 
Hessian matrix Hq^qi: 



where we have defined iJ? by 



^ 1 [ndTn] 

~ 2! sntsn r ' 



(2.16) 



(For later convenience, we have chosen a definition that 
differs from the standard one by a factor of 1/2.) More 
explicitly, we have 



(2.17) 



As it is adequate to be concerned with matrix elements 
to leading (i.e. first) order in e, we shall neglect higher 
orders from now on. 



Consider any variation {<5ri^.} of the {f^j.} around a 
stationary point {^j,} ■ To first order in {Jilj.}, the vari- 
ation of the fr ee e nergy functional is, of course, zero. We 
see from Eq. (2^) that the second order variation is |l^] 



^kik2k3 '^1 



, , A ft; Si}; 5Q; 

+k2 + k3,0 ki k2 ks 



(2.12) 



Now consider expanding around the liquid state for any 
value of e. In this case, the second variation reduces to 




(2.13) 



which evidently indicates that the liquid is stable for 
e < and unstable for e > 0. For e > the only can- 
didate we know of for a stable thermodynamic state is 
the amorphous solid. From now on we focus only on that 
state. 

By in serting the value of the order para meter , 
Eq. (2^), into the three- wavevector sum in Eq. ( 2.12 ), 
we obtain 



C. Change of basis 

In order to simplify the diagonalization of the Hessian, 
we are going to exploit the symmetries of the problem. As 
a direct consequence of the invariance of the free-energy 
functional under translations and the MTI property of 
the amorphous solid state, the matrix element Hj^ ^ of 

the Hessian only connects wavevectors k and / such that 
k = 1. This already reduces the complexity of the prob- 
lem by making the Hessian block diagonal. 

As Hf. I depends on A;^, r, and k ■ I, one might expect 
to find a symmetry under arbitrary rotations in {n + l)d 
dimensions, which would simplify the diagonalization of 
the Hessian still further. However, the factor is not 
invariant under some of those rotations. Instead, H only 
displays a rotational symmetry in nd dimensions, but this 
will enable us to simplify the task in much the same way 
as is commonly done for central potentials in quantum 
mechanics. 

In order to make this symmetry explicit, we choose a 
fixed matrix T G 50((1 -I- n)d) such that for any vector v 
in replicated space we explicitly isolate v from the other 
nd independent coordinates, which we call v: 



2 + fe3,0 "ki^O 3 ^ 



Tv : 



n+n 

V 



(2.18) 
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Due to T being orthogonal, scalar products remain sim- 
ple in the new coordinates: 



hrs 



V • w „ „ 

V ■ W = h V ■ w. 

l + n 



In the new coordinates, we have 



1 ^ k2 



2^ 1 



n 



(2.19) 



2e j d6'7r(6l)e-('^-')'/^^}. (2.20) 



This expression, taken naively, would immediately tell 
us that the Hessian is invariant under rotations in nd 
dimensions: 



GO 
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Irs 



H 

(original 
Hessian) 



hrs 




Vi? e 50 (nd) : 



(2.21) 



However, there is an important caveat. Our Hessian is 
only defined for wavevectors in the higher replica sec- 
tor, but the proposed rotations can take a vector in the 
higher replica sector and transform it into a vector in 
the one replica sector. For the moment we are going to 
ignore this difficulty, and simply diagonalize the matrix 



obtained by using Eq. (2.17) as its definition, with k and 
I taking any nonzero values, both in the higher and in the 
one replica sector. (This enlarged Hessian matrix will be 
termed the "extended Hessian" to distinguish it from the 
"original Hessian" which does not have those additional 
matrix elements. See Fig. |^.) After having diagonalized 
the extended Hessian, we will return to the issue of the 
one replica sector. For the moment, let us just anticipate 
that, in the replica limit n ^ 0, the only effect of this ex- 
tension of the Hessian on its spectrum of eigenvalues will 
be the addition of one spurious eigenvalue, corresponding 
to a fluctuation localized in the one replica sector. 



FIG. 1. Comparison of the original Hessian matrix (upper 
figure) and the extended Hessian matrix (lower figure). The 
extended Hessian matrix is depicted as being formed by four 
blocks, respectively connecting: the Irs with itself (labeled 
the Irs with the hrs (labeled H^"^), the hrs with the Irs 
(labeled H^^), and the hrs with the hrs (labeled H^^^). The 
block H^^ corresponds exactly to the original Hessian. 

For the modified problem of diagonalizing the extended 
Hessian, the SO{nd) symmetry holds, and Eq. (2.21) is 
correct without caveat. In what follows, we will use the 
same strategies as in the diagonalization of a quantum- 
mechanical hamiltonian for a particle in a central poten- 
tial. The role of the hamiltonian will be played by iJ^ ^. 
We will also exploit the symmetries of the problem: the 
rotational symmetry in nd dimensional space will allow 
us to write each eigenfunction as a product of a radial 
part, which will be obtained by solving a one dimensional 
eigenvalue equation, and an angular part, which will sim- 
ply be a surface harmonic function p9| ] in nd dimensions. 
The quantity k, which is exactly conserved by Hf, will 
play the role of a conserved quantum number. We will 
obtain distinct sets of eigenfunctions for each fixed value 
of k. 

We are going to work in the Hilbert space of complex 
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functions of the variable k (7^ 0). We exclude the origin 
because we are interested in fluctuations of the order pa- 
rameter, and JIq cannot fluctuate. We define the scalar 
product in this space as follows: 



(2.22) 



which simplifies, in the thermodynamic limit, to |21[] 



dk 



dk dk 



V 



r(k,fc)ff(k,A). (2.23) 



[l + nyi^ (27r)(i+")'' 
We define a basis set for this Hilbert space by 

^ppa (k, fc) — 

{l + nr/^{2nr'^/H-^^5m-p)p'^SAh)^ (2-24) 

Here {Sa{4>)} are the normalized surface harmonic func- 
tions defined on the unit sphere in nd-dimensional 
space jf^ , p is any wavevector in d dimensions, and 
p is any positive number. (Surface harmonics are ho- 
mogeneous trigonometric polynomials, they are general- 
izations of spherical harmonics to any space dimension 
d > 3. The label cr is a set of integers that charac- 
terize the appropriate trigonometric polynomial: for ex- 
ample, for d = 3, the surface harmonics are the usual 
spherical harmonics, and a = (l^m), with I the degree of 
the trigonometric polynomial and m a label that distin- 
guishes between polynomials of the same degree.) The 
notation ip^ = k/\k\ denotes the unit nd-dimensional vec- 
tor along the direction of k. The elements of the basis set 
Wppa} are o rthog onal and normalized under the scalar 
product Eq. ( |2.22| ): 



(<y5p'pV I'^ppcr) = 5{p' - p) (5p',p 5a' 



(2.25) 



As suggested above, we propose to express each eigen- 
function for the problem in the form: 

-(/'rp<T(k, k) 

dpRr{p)Vppa{kk). (2.26) 



where the discrete 5 function ensures that the eigenfunc- 
tion is localized on points with a fixed value of k, the 
surface harmonic gives the angular dependence on fc, 
and the radial function gives the radial dependence 
on k. By using the normalization condition for the basis 
set, we obtain the normalization condition for the radial 
part: 



dfc|i?^(fc)|^ 



(2.27) 



We now compute the matrix elements of the Hessian be- 
tween the elements of the basis set {(pppo-}- 

Consider first the part iJ^ of the Hessian that is di- 
agonal in fc, i.e., the first term on the r.h.s.of Eq. (2.20) 



This part is also diagonal in the basis {i^ppo-}, with the 
matrix elements 



{^p'f>'a'\H^\ipppa-) 

S{p' -p)Sp,^p Sa'^a 



(2.28) 



The non-diagonal part H'^ of the H essian, given by the 
second term on the r.h.s.of Eq. ( 2.2Cl| ), has, in the new ba- 
sis, matrix elements that only connect different values of 
the radial coordinate p but are still diagonal in p and a. 
It is shown in App. ^ that the matrix elements of H'-' 
have the form 



p'CT' 

<^p',p Sa 



iH^'liPppa) 



(-2e) C„ v\:Jip'/Ve,p/Ve), (2.29) 



with 



C„ = (e/47r)"'^/2(^^^)-d/2 



(2.30) 



and 

(n) /I \ r, / ; 

rji (x ,x) = zVx'x 

00 

31-nd/2 



f2x'x 

h-l+nd/2 I Q 



(2.31) 



where /^(x) is the modified Bessel function of order v. 
The label /(= \(t\) indicates the degree of the surface har- 
monic 5*0- as a trigonometric polynomial. The constant 
Cn satisfies the condition 



lim Cn 



1, 



(2.32) 



via which it disappears from the eigenvalue equation in 
the replica limit. The kernel rij^^ix' ,x) is real and sym- 
metric, and controls the non-diagonal nature of the ma- 
trix elements. Due to the positivity of Iv{y) for > — 1 

and 2/ > 0, the kernel vi["'\x' ,x) is positive for xx' > 0. 

For xx' 0+, vi["'\x' ,x) vanishes, except if Z = and 
nd> Q, in which case it is divergent. 

As both H'^ and are diagonal on the y> and a 
labels, the eigenvalue equation for the Hessian, 



can now be simplified to a radial equation: 



(2.33) 
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kR{p) 



dp' {(Pppa\H\Lpp'pa)R(p'), 



1^ 



Rip) 



1 + ny 

-2eCn J^"° ^^vl:l{p/Ve,P/Ve) R{p)- (2.34) 

This radial equation can be simplified further by the 
rescaling 



c 



p2 



p/V^, 



2(1+7 

u{x) e^^^Ri^x), 



(2.35a) 
(2.35b) 



which removes the e dependence from the eigenvalue 
equation, thus making both the eigenvalue C and the 
eigenfunction u{x) e- independent and obeying 



(u{x) ^ ^u{x) ~2Cn dx'r]l^hx,x') u{x'). (2.36) 
2 Jo 

In the replica limit, n 0, this radial equation reduces 
to 



(^u{x) = ^uix) - 2 / dx' ril'^hx,x') u{x'). (2.37) 
2 Jo 



For all cases except jfjl = this limit is straight- 
forward, because r]'^^^^ (x,x') smoothly converges to 
■q^^'l {x , x' ) . For the special case of \a\ = 0, the limit 

n for r]^^^^ (x, x') is singular near the origin. Let us 
mention here a property of Eq. (2.37) that does not ap- 
ply to Eq. ( ]2.36D . As /-_i(z) = /i(z) for all values of the 
variable z, we have the equality 



Vo ) = 77^ '{x,x ) 



(0), 



(2.38) 



which means that the radial equation (2.37) is the same 
for \a\ = and for \a\ ~ 2. In Sec. IIIC we discuss 



in more detail the relations between the solutions to 



Eq. ( peD and Eq. (|2^). 

Both radial equations, Eq. ( 2.36| ) and Eq. (2.37), are 
eigenvalue equations for Hermitian operators. This guar- 
antees the existence of a complete orthonormal basis of 
eigenfunctions, all of them having real eigenvalues. No- 
tice also the nontrivial fact that the radial equation is 
well defined in the replica limit ri ^ 0. 

The form of the radial eigenvalue equation tells us that 
the radial eigenfunction and the eigenvalue depend on 
the degree I = \a\ of the surface harmonic considered 
(which plays a role analogous to that of the total angular 
momentum quantum number / in the central potential 
problem for a quantum mechanical particle), and on an 
additional label r, playing a role analogous to the ra- 
dial quantum number in quantum mechanics. Therefore, 



the eigenvalues of the extended Hessian are given by the 
relation: 



Kir{k) = (1 + Clr)e + 



(2.39) 



As it is easier to work with scaled variables, let us ex- 
press the condition that there be no unstable fluctuation 
directions (i.e., k > 0) in terms of C- 



K;r(k) > 



Vk 



Or + 1 > 0. 



(2.40) 



The right hand side of the equivalence sign is the condi- 
tion that we are going to establish in what follows. 



III. LANDAU THEORY: EIGENVALUES OF THE 
HESSIAN 

In this section we are going to establish positive lower 
bounds for all the eigenvalues of the Hessian of the Lan- 
dau theory, except for the zero mode associated with the 
spontaneously broken translational symmetry. In other 
words, we study the set So containing the limits, when 
n —^ 0, of the eigenvalues of the original Hessian. How- 
ever, for technical reasons, it is convenient to first study 
two other sets of numbers, denoted by Se and Sr- Se is 
the set containing the n — > limits of the eigenvalues of 
the extended Hessian; eac h elem ent in Se can be written 
in the form given by Eq. ( 2.39| ), where Cir is taken to be 
the n ^ limit of an eigenvalue in Eq. ( 2.36| ) . Sr is the 
set containing all numbers k;,- (k) computed according to 
Eq. (2.39), wit h Cir chosen to be an eigenvalue in the 
radial equation (2.37). 

In Sec. Ill A we show that Sr contains a zero element 



corresponding to the zero mode associated with the spon- 
taneously broken translational symmetry. In Sec. HI B 
we compute positive lower bounds for all other elements 
ofSr- In Sec. IIIC, we show that the only difference be- 



twee n the limit for tt. ^ of the eigenvalue spec trum of 
Eq. ( ^.36| ) and the eigenvalue spectrum of Eq. ( 2.37 ) is 
that in the former a spurious eigenvalue corresponding to 
fluctuations in the Irs appears, which is not present in the 
latter. Therefore the spectrum Se of the extended Hes- 
sian contains a spurious eigenvalue not present in Sr- We 
also show, in Sec. HI C , that the eigenvectors of the orig- 



inal Hessian correspond to all the eigenvectors of the ex- 
tended Hessian, except the spurious one, i.e. that S p and 
Sr a re ide ntical. Using the results obtained in Sees. [II A 
and IIIB| , this will allow us to conclude that the amor- 
phous solid state is locally stable. 



A. Obtaining the zero mode 

We first consider the eigenfluctuation associated with 
the translational symmetry, and show that it is a zero 
mode. From Eq. ( p.9[ ), we see that this fluctuation can 
be written as 



7 



(3.1) 



The only angular dependence of (5f2^, is given by the pref- 
actor k-d, which is a degree-one polynomial in k. This 
guarantees that this fluctuation resides in the |(7| = 1 sec- 
tor. By taking the scalar product with the appropriate 
element in the basis {v'ppcr} [which we label by ct = (1, 0) 
by analogy to the spherical harmonic Yio oc z/r], we ob- 
tain the radial function associated with S^f.'. 



R{k) = ((/3p,p=o,fT=(i,o)l'5^^>, 
Jo 



(3.2) 



where An is a numerical prefactor, which we can ignore in 
what follows. Taking the replica limit, and transforming 
to scaled variables, we obtain the scaled radial function 



u{x) 



IX I d9T:{9)e 
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(3.3) 



In App. (^) we show by explicit computation that this 
form for u{x) satisfies Eq. ( 2.37 ) with C, = —1. By 
Eq. (2.39p, this means that the corresponding 5U,f^ given 
by Eq. (3.1) is an eigenvector of the Hessian with zero 
eigenvalue. [As SVLf, given by Eq. (3.1) is only nonzero 
for k in the hrs because otherwise fc = 0, it is simultane- 
ously an eigenvector of the extended Hessian and of the 
original Hessian.] 



B. Positive lower bounds for the eigenvalues 

Having obtained the zero mode for a specific form of 
fluctuation, we now discuss generic fluctuations in the 
order parameter field, and show that all the other eigen- 
values are positive-definite. For the case ^ 7^ 1, we will 
obtain positive lower bounds for t he eigenvalues by an- 
alytical manipulat ion o f Eq. ( 2.37 ). For the case / = 1, 
we will solve Eq. (2.37) numerically and show explicitly 
that the lowest eigenvalue corresponds to the zero mode 
already obtained, and that all other eigenvalues corre- 
spond to positive values of K(k). 

Consider one p articular scaled radial eigenfunction 
u{x) in Eq. ( |2.37| ), with eigenvalue C,. To simplify the 
argument we temporarily switch to the normalization 



dx \u{x)\ = 1, 

and we define the quantity 

s{x) = sgn[u{x)). 



(3.4) 



(3.5) 



In what follows, we express the eigenvalue C in an un- 
usual but convenient form that allows a lower bound to 



be de rived from it. By combining the eigenvalue equa- 
tion (2.37), with the normalization condition (3.4) and 
the definition (|3.5|), we obtain 



( — C dx \u{x) 



dx s{x) C u{x) 



dx s{x) —u{x) — 2l dx' ri^^{x,x')u{x') 



dx—\u{x)\ 



-2j dxdx' r]^°^^{x',x)\u{x)\s{x)s{x'). (3.6) 

(In the last line we have interchanged the dummy vari- 
ables x and x' .) From the expression just derived for 
the eigenvalue Ci it follows from the nonnegativity of 
Vi°]ix',x) that 



C > / dx-fi„i{x) \u{x)\ > 7|^|, 
Jo 



X 



7/(x) = y-2y^ dx'^f\x',x), 
7/ = inf 7/(a;). 



(3.7a) 

(3.7b) 
(3.7c) 



Here, the symbol inf indicates the greatest lower 
bound [^ for a set of real numbers. 

It is convenient to write 7; (a;) in terms of another func- 
tion /3i(u), as follows: 



(3.8a) 



^i{x) = / deTT{e)i3i{x/Ve), 

Jo 

-4 / rf^i^^^e-("'+^')/^_l(2^t^;). (3.8b) 



Here, we have used the definition of an average with re- 
spect to the distribution tt. 



de7T{0)f{e), 



(3.9) 



and we need, in particular, the numerical value 

{0-'^)^ w 0.881768, (3.10) 

which can be obtained by using the function 7r(0) of 
Ref. As tt{6) is nonnegative and normalized to unity, 
Eq. (p.8a[) implies that 



C > 7|<t| > I3\a\, 



(3.11) 



where 
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A =inf/3/(z;). 



(3.12) 



5.0 



The bounds for different values of the index / are not 
independent. In fact, because Ii{x) < Ii-i{x) for a; > 
and ^ > 1, we have the inequahty T?;°\(a;, a;') < 77;'-*''' {x,x') 
for xx' > and I > 1, and from this inequahty it follows 
that all the bounds defined so far [jiXx), 7/, Pi{v), and 
/?/] are increasing functions of I for I > 1. Thus, if we 
obtain a positive lower bound for one value of ^ > 1, the 
same bound applies for all larger values of I. 

In order to obtain more concrete results, we need an 
explicit expression for I3i{v). In App. (|^) we obtain the 
exact expression 



r( 



(3.13a) 



along with the asymptotic forms j3'^^\v) and /3^^-'(u), 
given by 



for w > 1, (3.13b) 



for w < 1, (3.13c) 



as well as the lower bounds 

I3i{v) > I3^>\v) 



for / > 1, 



(3.14a) 
(3.14b) 



1 = 1 




-1.0 



0.0 



3.0 



1.0 2.0 

V 

FIG. 2. Plot oil + l3i{v) (for < Z < 4) and 1 + /3<>' («) as 
functions of v 

Let us now obtain the lower bounds 1 + (3i for 1 + Ci 
and show that they are positive for ^ 7^ 1. For I = 4 (and, 
as /?/(w) grows with I, for all ^ > 4), /5/(w) is positive for 
all nonzero u, and thus Pi = /3;(0) — 0. For < / < 3, 
[3i is obtained by numerically minimizing Eq. (3.13a). In 
Table ||, we give the numerical values for these bounds. 
These lower bounds establish that all of the K;r(k) are 
positive for 1^1. 



TABLE L Lower bounds for eij 


jenvalues of the Hessian. 


/ 


1 + A 


1 


-0.55571 


0,2 


0.27376 


3 


0.94274 


> 4 


1 







Here T{z) is the Gamma function, and M(a, 6, z) is a 
confluent hypergeometric function ( |^^, Chap. 13) 

As mentioned above, we need to show that 1 + C is pos- 
itive. Thus, the quantity of interest is really 1 + f3i{v), as 
opposed to l3i{v). In Fig. || we plot 1 + Pi{v) as a func- 
tion of V for < ^ < 4, together with its asymptotic form 
1 + /3^^^(w) valid for large values of the argument v. 
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Let us now focus on the remaining sector, namely 1 = 1. 
As, for this case, our lower bound is negative, we can not 
draw any conclusion from it. We have already shown that 
there is a zero mode, but there could still be one or more 
negative eigenvalues, which would render the proposed 
amorphous solid state unsta ble. The numerical solution 
of the radial equation (2.37) for this case yields, within 
numerical error , the following two lowest eigenvalues 
(both nondegenerate): 



1 + Cio = -0.00002 ± 0.00005 
1 + Cii= 0.98412 ±0.00009 



(3.15) 



Evidently, 1 + Cio corresponds to the expected zero mode, 
and we can conclude that there are no further zero modes 
and that all other eigenvalues are positive definite. 

To summarize, we have shown that 1 + C ^ for any 
eigenvalue C of Eq. ( ^.37| ), i.e., that all elements of Sr 
are positive or zero. In the next subsection we will show 
that Sr and So are identical, and therefore that all of the 
eigenvalues of the original Hessian are either positive or 
zero. 



C. The one replica sector and the spurious 
eigenvalue 

We now need to return to the issues that we postponed 
earlier, namel y ou r extending of the Hessian matrix de- 
fined by Eq. (2.17) so that it also be defined in the one 
replica sector, a nd th e diffe rence s between the spectra of 
the radial Eqs. (|23|) and ( p3^ ). 

As the Hessian matrix (both in its "original" and its 
"extended" versions) leaves exactly uncoupled fluctua- 
tions with different values of k, it is consistent to con- 
sider separately the MTI fluctuations (those with k = 0) 
and the non-MTI fluctuations (those with k ^ 0). 

For the case of non-MTI fluctuations we will show that, 
in the limit n ^ 0, the hrs and Irs are not coupled by 
the extended Hessian matrix. Furthermore, each one of 
the eigenvectors belongs to one of the sectors, inasmuch 
as it has a vanishing overlap (as n — > 0) with vectors in 
the other sector. 

To understand this issue, we need to look at the form 
that the Irs and hrs take, in the replica limit. For a 
wavevector 



p= (0,...,0,p,0,...,0) 

in the one replica sector, we have 

p p and — P^, 



(3.16) 
(3.17) 



and, by using Eq. (2.19) with v = w = p, we have 



\P\ ^ \ P 



1 + n 



1 + n y 1 + n 



(3.18) 



This means that the radial coordinate \p\ goes to zero 
like n^/^ in the replica limit, n — > 0. Moreover, for each 
fixed value of p (= p), the n + 1 wavevectors defined by 



e„(p)^(eO,ei,...,e"), 

for /3 7^ a, 
p for = a, 



(3.19a) 
(3.19b) 



with a = 0, . . . ,n, are the only vectors in the one replica 
sector that satisfy the condition that the sum of their 
(n + l) component d-dimensional wavevectors is equal to 
p. These two results tell us that the one replica sector 
corresponds, for fixed p, to a set of n -I- 1 points that, 
in the replica limit, converge to the origin of p space. 
Consequently, to see whether or not a given eigenvector 
has any overlap with the Irs, one needs to look at the 
properties of the corresponding radial eigenfunction near 
the origin. 



Let us then consider the scaled radial equation (2.3(;) 
for the region close to the origin, and let us keep n > for 
the moment. By using the small- argument behavior of 
the modified Bessel function (vahd for ^ —1, —2, . . .), 



r(^ + i)' 



(3.20) 



we obtain the asymptotic form of the kernel 77 for x <C 1 
and x' 5; 1: 



miy) 



X 



l + (nd-l)/2 



r{l + nd/2) 

o 

d0Tr{e)0' 



mi (x'). 



I q-vVS yl+{nd-l)/2_ 



(3.21a) 
(3.21b) 



By inserting this asymptotic form into Eq. ( p. 36 ) , we ob- 
tain 



2-2 _42,'+(nd-l)/2 

{C~^)u{x)^ r(/-fnd/2) 



Ui 



(3.22a) 

dymi{v)u{y). (3.22b) 



For n positive and small, Eq. (3.22a) can only be satis- 
fied for C 7^ The term proportional to x'^ is thus 
negligible, and we obtain, for x ^ 1, 



u{x) 



T{l + nd/2) 



(3.23) 



The leading behavior of this radial eigenfunction for nd 
small and positive depends on the value of the degree 
I = \a\ of the surface harmonic function. For I = 
there is one eigenfunction that diverges at the origin like 
^(nd-i)/2 eigenvalue C- is given by the expression 
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which, in the rephca hmit, becomes 
hm C- = -2. 



(3.25) 



The presence of this divergent eigenfunction as a solu- 
tion of Eq. ( ^.36| ) depends crucially on the singularity of 

77g"^(a:, x') at the origin for n small but positive. It is not 
a solution of Eq. ( 2.37 ), which only has solutions that go 
to zero at the origin: by Eq. ( 2.3§| ), the radial equation 
for n = and Z = 0, is identical to the radi al eq uation 
for n — and 1 — 2, which means, by Eq. ( 3.23| ), that 
the radial eigenfunction satisfies the condition 



\uix)\ 



(3.26) 



We will show below that the divergent eigenfunction 
present for n > corresponds, in the replica limit, to 
an unphysical fluctuation, i.e., a fl uctu ation i n the one 
replica sector. In fact, from Eqs. (2.39) and (3.25), we 
see that its eigenvalue K_(k) is negative for small k: 



«_(k) 



-e + 



(3.27) 



However for I = 0, all other radial eigenfunctions are 
orthogonal to the one just found, and thus make the in- 
tegral Ui vanish. Their behavior is controlled by the next 
power in the expansion of 77q"^(x, x'), and consequently 
they vanish for x ^ at least as fast as a;("'i+3)/2^ More- 



over, for Z > 0, by Eq. (3.23) we see that all the radial 
eigenfunctions vanish for a; ^ as or faster. 

These results can be summarized as follows: for n 
0+ all but one of the radial eigenfunctions are regular at 
the origin. The one singular eigenfunction corresponds 
to / = and scales like x'"''"-'^'/^ for x ^ 1. The regu- 
lar eigenfunctions can have any value of / and vanish for 
X ^ as or faster. 

In all cases in which the eigenfunction is regular at 
the origin, it is permissible to take the limit n — > in 
Eq. (2.36). This is because these eigenfunctions van- 



ish a t the origin fast enough that the integral term in 
Eq. ( 2.36 ) does not pick up any extra contribution from 
the singularity of r]l^^\x,x') [which, by Eqs. ( 3.21a ) and 
( |3.21b| ), is at most of order nd{xx'Y'^'^-^^/'^]. Thus, the 
spectrum of eigenvalues of E g. (^.37 ) is the same as the 
limit of the spectrum of Eq. ( 2.36 ) when n ^ 0, except 
that the spurious eigenvalue C- is absent in the former 
and present in the latter. 

We now show that the one replica sector fluctuations 
decouple from the higher replica sector fluctuations in the 
replica limit. Consider the following complete orthonor- 
mal basis set for the fluctuations in the one replica sector 
with fixed p = p: 



w 



(k) = X!^J."'^fc,e„(p)' 
a=0 

Vn+l 



i2TTja/{n+l) 



(3.28) 



where j = 0, . . . , n. 

Let us compute the scalar product (wjliprpa) of one of 
these basis functions for the one replica sector with one 
of the eigenfunctions of the extended Hessian, which has 
the general form given in Eq. (2.26): 



{Wj\lprpa 



-^^u;;(^)^Arp.(fc) 



■{V^\P\) Rr{V^\P\) 



(3.29) 



Here, we have made use of the relation Eq. (3.18). There 
are two possible cases to consider, depending on whether 
or not Rr is singular at the origin. If it is singular, we 
have / = and, for small k, 



Rik) 



(3.30) 



wher e M is a normalization constant determined by 
Eq. ( 2.27 ). Its value is given by 



= Vnd{l + 0{nd)) . 



(3.31) 



As / = 0, the angular part of iprpa is isotropic, and is 
given by 



(3.32) 



where Tnd — 27r"'^/^/r(n(i/2) is the surface area of a 
unit s p here i n nd-di mcns ions. By combining Eqs. (3.29), 
(|3^, ( |33l| ), and (|33|), we obtain 



Q=0 



(5,-o(l + 0(«)) 



(3.33) 



This result implies that, in the limit n ^ 0, the eigen- 
function that is singular at the origin lies entirely in the 
one replica sector. 

Let us now consider the case in which Rr is not sin- 
gular at the origin. In this case, for small k, the radial 
eigenfunction has the form 



R{k) 



£-1/4 li(fcM) < AAe-(l'-i|+i)/2yfc(l'-i|+i/2) , 



(3.34) 



where the normalization constant Af does not vanish in 
the replica limit. As, in this regular case, I need not be 
zero, we have to obtain an estimate for the normalization 
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constant of the surface harmonic for all values of /. Con- 
sider a monomial Mm{<t>) defined on the ZJ-dimensional 
unit sphere 



MM = 4>T ■■■(!>[ 



nil 
mo _ ^1 



(3.35) 



Here (xi, . . . , xd) are the cartesian coordinates of a point 
X, r = {x\ + • • • + a;|,)^/^ is the radial coordinate for the 
same point, and (j) = {(pi, . . . , ^u) = x/r is the unit vec- 
tor pointing in the direction of x. The integral of the 
monomial over the unit sphere is 



.M™((/.) = 



(fx X 



mi 
1 



■ x^J^" e^^'^i^ — ^^^^ 



^rniA [-mo 



2nr(i 



= < 



if rrij even Vj , 



otherwise. 



(3.36) 



In the case of interest to us D ~ nd and X^^i "b' = 2|(t| 



From Eq. (3.36) we conclude that the normalization fac- 
tor No- for the surface harmonic So- has, in the n — *■ 
limit, the asymptotic form 




1/2 for \(7\ = 0, 
for \(7\ ^ 0. 



(3.37) 



Here, we have ignored factors that have finite limits when 
n 0. This result can be summarized as follows 



via) 



-1/2 for \a\ = 0, 
for |ct| ^ 0. 



(3.38) 



By inserting Eqs. (|3.34|) and ( |3.38|) into Eq. (|3.29|) , we 
obtain the following scaling with n for the sought scalar 
product: 



1 " 



w w 

2 J m,Q :?,CK 



-2e / dOnid) e-[^°(P)-^''(P)l'/^^<,,^,-^' 



^2 , 



,+ E_)+0(n). 



(3.40) 



Thus we see that, as expected, the eigenva lue o btained 
here is the same as the one obtained in Eq. (3.27) for the 
singular eigenfunction of the extended Hessian. 

In summary, for non-MTI fluctuations, in the replica 
limit all regular eigenfunctions of the extended Hessian 
are orthogonal to all of the Irs vectors, and the singular 
eigenfunction of the extended Hessian coincides with the 
isotropic (j = 0) fluctuation in the Irs. Consequently, in 
the replica limit, the higher replica sector is an invariant 
subspace for the extended Hessian, and therefore the reg- 
ular eigenfunctions of the extended Hessian are the eigen- 
functions of the original Hessian. More significantly, the 
eigenvalues of the original Hessian are the eigenvalues of 
the extended Hessian for its regular eigenfunctions. 

For the case of MTI fluctuations, their components in 
the one replica sector are exactly zero, because the con- 
ditions k G Irs and k = are incompatible. However, all 
of the argument presented above still holds, except that 
now the radial eigenfunction that is singular at the origin 
coincides with a spurious fluctuation in the zero replica 
sector (i.e., a fluctuation of Jig) The spectrum of 

the original Hessian is, also in this case, given by the 
eigenvalues corresponding to radial eigenfunctions regu- 
lar at the origin. Thus, we have shown that the sets 
Sr and S„ are identical, and we can use the results of 



Sees. [II A and IIIB to characterize the spectrum of the 



original Hessian. 

The eigenvalues of the original Hessian have the gen- 
eral form 



2 ' 



(3.41) 



with 



\{w,\^rp.)\ < (yH|k|)V2(VH|k|)l/2+l'-l| n-M, 



< 



.1/2 



for / = 0, 



|i-l| + l)/2 



(3.39) 



This relation shows that in the limit n those radial 
eigenfunctions that are regular at the origin give rise to 
eigenvectors that lie entirely in the higher replica sector. 

For completeness, we now also compute explicitly the 
matrix elements of the extended Hessian between mem- 
bers of the basis set {wjYj^Q for the one replica sector 
fluctuations with k = p. 



1 + Cio = 



and 



1 + C/r > 



(3.42) 



for (Z,r) ^ (1,0). (3.43) 



Therefore there is a zero mode corresponding to [l, r) — 
(1,0) and k = 0, which is continued by a branch of soft 
modes with eigenvalues 



Kio(k) = — (> for k ^ 0). 



(3.44) 
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All other eigenvalues are positive, with one continuous 
branch of modes labeled by k for each value of {l,r). 
The minimum eigenvalue for each branch is given by 



KiriO) = (l + Or.)e > 0, 



(3.45) 



which goes to zero as the transition is approached (i.e. as 
e 0). Consequently, the amorphous solid state of 
Ref. is locally stable near the transition. 



Landau theo ry dis cussed above, i.e., it is also described 
by Eqs. (H), dU), and 

We now expand the free energy functional to quadratic 
order around a stationary point, and obtain its second 
derivatives with respect to the fields {i^g}. In this section 
we use the notations H and H to refer to the exact Hes- 
sian for the microscopic theory and the extended Hessian 
for the Landau theory, respectively. For k and k' both in 
the higher replica sector we have 



IV. RANDOMLY CROSSLINKED 
MACROMOLECULES 

We now consider one example of a semi-microscopic 
theory that exhibits the amorphous solidification tran- 
sition, namely the case of randomly crosslinked linear 
macromolecules. In this theory there appears a control 
parameter /i^ (= 1 + e/3) that determines the crosslink 
density, and such that the system exhibits the liquid 
phase for < 1 and the amorphous solid phase for /i^ > 
1 . In this semi-microscopic theory, the field flp (with p in 
the one replica sector) is present and allowed to fluctuate, 
and there is a couphng parameter = — fM^ jj- asso- 
ciated with its fluctuations. (The parameter gives the 
strength of the excluded-volume interaction between the 
macromolecules.) The free-energy functional (per macro- 
molecule) has the form ||l|,Q 



- In /exp (iK^Yl -^^ ^pPi 



V 

yn kf'h 



W 



n+1 



(4.1) 



where the symbol denotes a sum over replicated 
wave vectors in the one replica sector, and the f symbol 
additionally restricts any summation to the half space 
of relevant wave vectors [i.e., d-dimensional or (n-l-l)d- 
dimensional] such that their scalar product with a fixed 
unit vector (n or n) is positive. Here, we have used the 
definition of the one-macromoleculc Fourier transformed 
density Pj. , i.e.. 



Pf. = / ds cxpik-c{s), 
Jo 



(4.2) 



for a macromolecular configuration c(s), and the Wiener 



replicated average 



^ is defined by 



(O) 



/PcOexp{-i/>|^f} 



ri+l 



JVC exp{-i/>|^| } 



(4.3) 



Let us note here that to leading order in e the amorphous 
solid stationary point in this theory is the same as in the 



k,k' 



2 / 2 
yn \ "kM ^irAP-kPk'ln+1,1- 



k'^\ 



deTr{9)e- 



'(k-k'f/ee 



3 kl' 



(4.4) 



For k in the higher replica sector and p in the one replica 
sector we have 

_ ■\2^[}^/ \W,n 



■ yi+n 

rxiNp' 



V 



3 



d0Tr{0)e- 



-^n-"^ ^hl 



(4.5) 



Finally, for both p and p' in the one replica sector we 
have 



6'^[ndTn] 



= H, 



p,p' 



V V P'^' 

Op.p' y \ 



XI N 



1 



V 
V 



\P-pPp')n+l.c 



[l + 0{e) + 0{f)]]. (4.6) 



In obtaining these formulas we have made use of the def- 
inition 



( O exp ( 



w 

n+l 



exp(^A2fEp^^pP; + '(^E^%.P^ 



w 

n+l 

(4.7) 



The notations iJ'^'^, H^^ , and -ff^^, respectively, refer to 
the higher replica, cross-sector, and one replica parts of 
the Hessian matrix. Fig. |3| depicts the relation between 
H and H. 
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hrs 
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FIG. 3. The Hessian matrix H for the semimicroscopic 
theory in terms of the extended Hessian matrix H for the 
Landau theory. Each block of H is written in terms of the 
corresponding block of H, except for H^^ . {H^^ is nonzero 
when e — > 0, as opposed to H^^, which vanishes linearly with 
e near the transition.) 

As MTI (i.e. k = 0) fluctuations do not have any com- 
ponent in the Irs, the relevant Hessian in this case is just 
i^hii = (1/3)5^'^. Consequently the results obtained for 
the Landau theory tell us that there is a zero eigenvalue 
corresponding to the anticipated zero mode, and that the 
remaining eigenvalues are positive. 

Let us now consider general fluctuations. We will show 
that, in the replica limit, the eigenvectors of the Hessian 
H for this problem are the same as the eigenvectors of 
the extended Hessian H for the Landau theory, and the 
one replica and higher replica sectors are again invariant 
subspaces for this Hessian. 

Let us consider a regular eigenvector \iprpa) of H and 
one of the elements of the basis set ||w ,-)}"^n of the one 



replica sector fluctuations. By Eq. (4.6) 



H 



11 1 



''Irs 



V 



1 



V ' 



Oie) + 0{f)]}, (4.8) 



and therefore 



\{w,\H''\yJrf.a)\ = |«lrsKl^rp.)| < O(V^). (4.9) 



Analogously, by Eq. ( |3.40| ), 



H''\w,) = {- e+^]\w,)+Oin), (4.10) 



{w,\H''\4,rpa) 



(4.11) 



By combining Eqs. (4.5), (4.9), and (4.11), we can now 
estimate the matrix element: 



^(z«,|il"' + ff"|^.p.)+0(V^) 



< 



(4.12) 



This means that, in the replica limit, H\iprpa) has no 
projection in the one replica sector, and also that H\wj) 
has no projection in the higher replica sector. Therefore, 
also in this problem the one replica sector and the higher 
replica sector are decoupled invariant subspaces of the 
Hessian in the n — > limit. In the one replica sector, 
the eigenvalue is k^j.^ > 0. In the higher replica sec- 
tor, as H^^^ — the eigenvectors are the same 
as for the Landau theory, and the eigenvalues are ob- 
tained from those in the Landau theory by multiplying by 
1/3. As discussed before, all of these eigenvalues are pos- 
itive, except for a unique zero mode. Thus, also for the 
semi-microscopic theory of randomly crosslinked macro- 
molecules, the amorphous solid state of Ref. ||] is locally 
stable near the transition. 



V. SUMMARY AND CONCLUDING REMARKS 

In this Paper we have shown that in a system with ran- 
dom constraints near the liquid-amorphous-solid transi- 
tion, the amorphous solid state of Ref. is a locally 



18|. In order to do this, 
ue spectra of the stability 



and 



stable thermodynamic state [ 
we have examined the eigenva 
matrices, in the contexts of both the Landau theory for 
the transition and a semi-microscopic model of randomly 
crosslinked macromolecular systems. In both cases the 
spectrum turns out to be non-negative, with only a sin- 
gle zero eigenvalue, and all the others positive. 

Let us remark that even though we do find a zero eigen- 
value for the stability matrix, we still declare that the 
stationary point is locally stable, as opposed to locally 
marginally stable. This is because in this system transla- 
tional invariance is spontaneously broken, and therefore 
there is a manifold of equivalent states that have exactly 
the same free energy and are connected to each other by 
the continuous symmetries of the system. The zero eigen- 
value (a.k.a. Goldstone mode) simply indicates that the 
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free energy does not change if one applies an infinitesimal 
translation to the thermodynamic state. 

In close analogy to the phonon spectra of ordinary 
solids, the fluctuation eigenvalues can be classified into 
two types: a soft branch of modes associated with "al- 
most rigid" displacements of the whole system (analo- 
gous to the acoustic phonon branch), with eigenvalues 
Kio(k) — and a set of stiff modes in which the 

structure of the system is altered more strongly (analo- 
gous to the set of optical phonon branches), with eigen- 
values Kir(k) = e(l -I- Cir) + In addition, there is in 
our case a softening of the system, because the eigenval- 
ues of the stiff modes go to zero at the transition. 

We have only addressed the issue of the local stability 
of the amorphous solid state. It is much harder to deter- 
mine whether the amorphous solid state is globally sta- 
ble, as the order parameter space to be explored is enor- 
mous. In particular, one could consider the possibility 
of a replica-symmetry breaking saddle point also being 
present and dominating the physical behavior of the sys- 
tem However, there are strong indications (although 
by no means conclusive evidence) that the (replica sym- 
metric) saddle point considered here is indeed globally 
stable. These indications mainly come from molecular 
dynamics simulations as the solid state observed in 
the simulations appears to be identical to the one pro- 
posed in Ref . ^ . 

An intriguing problem, left open for further study, is 
to establish how the structure of the eigenvalue spectrum 
of the Hessian matrix, and in particular the softening of 
the system near the transition, manifest themselves in 
the dynamics of the system. 
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APPENDIX A: NON-DIAGONAL MATRIX 
ELEMENTS FOR THE HESSIAN 

In this Appendix we collect some useful information 
concerning surface harmonic functions, and use it to com- 
pute the matrix elements of the non-diagonal part H'^ of 
the Hessian in the basis {^ppa}- 

The Gegenbauer (also called hyperspherical) polyno- 
mials play a role in D dimensions and with regard to 
the surface harmonics So- analogous to the role Legendre 
polynomials play in 3 dimensions and with regard to the 
spherical harmonics Yim- The Gegenbauer polynomial 
Ci{x) of degree / is defined by the generating function 
(see, e.g., Ref. |o|, Vol. II, Sec. 11.1.2) 



(1 - 2a;i + i^)-" = ^ CHx) tK 



(Al) 



i=0 



There is a generalization to dimension = p + 2 of the 
addition theorem for spherical harmonics, which relates 
the Gegenbauer polynomial to a sum of surface harmon- 
ics (see, e.g., Ref. Sec. 11.4): 



\a\=l 



h{l,p) 
4^1+p/2 



(2« + p)r(p/2) 



Y,S*M')S.{4>). (A2) 



\a\=l 



Here, (j)' and are any unit Z?-dimensional vectors, \a\ 
is the degree of the surface harmonic Sa as a trigono- 
metric polynomial, h{l,p) is the number of linearly in- 
dependent surface harmonics of degree I in dimension 
p -f 2, and tb = 27r^/-^/r(£'/2) is the surface area of 

1 /2 

a Z3-dimensional unit sphere. [As C; (x ) is equal to 
the Legendre polynomial Pi{x), formula (A2) reduces, 
for D = 3, to the usual addition theorem.] 

We also make use of the identity (see, e.g., Ref. [ pO[ , 
Vol. n, Sec. 7.15) 



2"^r(z.)^(n + i^)C;:(x)/„+.(z), 



(A3) 



n=0 



where /,y(z) is the modified Bessel function of order ly. 

In the case of dimension D — nd, by combining 
Eqs. (AS) and (A3), the following identity is obtained 



exp(x0'-0) = 27r"''/^(a;/2)i-"'^/2 

OO 

J2s*M)SA<t>)- (A4) 

/=0 \a\=l 

Here, x is any real number, and 0' and are unit nd- 
dimcnsional vectors. 

Let us now compute the matrix elements of the non- 
diagonal pa rt of th e He ssian in the basis {ifppa}- By 
using Eqs. dlH) and (^!20|) we obtain 



{(Pp'p'a'\H°\(Ppp^) 

, o f dkdldkdl nj/9/„ n,,^ , i-"*^ 

y (l-fn)^(2.)a.-)^^ ^^+-^ ^ ^^-^^ ' 

x5-,£<5(|fc| -p')^:'(0£) p'^S.jSi\l\~p)S,i^j) 



X %i (-2e) / (i6l7r(6i)e-(^-')'/^ 



(A5) 



-2e (27r)-" 



k,l 



k.l 



X / 6{\k\^p')\kr''-'d\k\ / s{\i\~p)\ir-'d\i\ 
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rOO 

Jo 

X / dcf>^ I d(t>^ SI, eM'^PP'h-H/^^) S.iH)- (A6) 



In the second step, we have separated the k and / inte- 
grals into radial and angular parts. The angular inte grals 
can be performed with the help of the identity (A4) and 
by using the orthonormality of the surface harmonics to 
obtain 



(_2e)eM-i)/2 



(^p'-^+p^y,e r ('ip'p\ 

(A7) 



which is equivalent to Eqs. ( |2.29D , (U^) and ( ^.3l[ ) 



APPENDIX B: RADIAL EQUATION FOR THE 
ZERO MODE 



dx' rif^x^x' ) u[x') 



-2Vx / dOde 
Jo 



/ d^e^"'/^ / d0'TT{0)Tr{e') 



(B4) 



Finally, we combine Eqs. (Bl) and (B4) to obtain 
u{x) + —u{x) — 2 J dx'rj^^'l^x^x') u{x') 

+ {l~e)TT{9)- j d^'7r(6i)7r(6'')} 

= 0. (B5) 

The justification of the last equality comes from the fac- 
tor i n b races being zero by the stationarity condition, 
Eq. dJI). 



Thus we have shown that Eq. (2.37) is satisfied by u(x) 
with the eigenvalue ^ = — 1. 



In this A ppen dix we show that the scaled radial func- 
tion of Eq. (3.3) corresponding to a change in the system 
due to a rigid displacement is a solut ion of the scaled 
radial eigenfunction equation Eq. (2.37) with — ~\. 

Let us first consider the diagonal term. By inserting 
the explicit form for u{x), and then performing an inte- 
gration by parts we obtain 



dOe- 





oo 



de 



roc n2 J 

-'V^l dee~^ /%--n{e) + en{e)}. (bi) 

Now the non-diagonal term gives 

/•oo 

-2 / dx'Tjf'\x,x') u{x') 
Jo 

= -2 rdx'd9d0'2V^^e-^^'^^"y^IoC-^) 
Jo ^ \ J 

By making use of the identity |^9[] 



(B2) 



2 1 dxxe 
'o 



loibx) 



(B3) 



we perform the integration over x' in Eq. (B2), thus ob- 
taining 



APPENDIX C: COMPUTATION OF LOWER 
BOUNDS 

In this Appendix we study in detail the bound-function 
Pi{v). We decompose Pi{v) as follows: 



2(0-1) 



ji{v) = 2 du^/^e'^'' hi^i(2uv) (CI) 

^0 

We now compute analytically the integral defining ji [v) : 

3i{v) = V^e'^" rdyy~^'^e~yii^^{2v^) 



r(V2 + i/4)^._,,,^...^ 



r(0 



2 4' 



where AI(a, b, z) is a confluent hypergeometric func- 
tion ( ||2^ , Chap. 13). B y inse rting this expression into 
Eq. (^), we obtain Eq. (|l3|). 

We can obtain more information by using the following 
integral formula for the confluent hypergeometric func- 
tion ( m, Chap. 13), valid for Rea > and Re6 > 0: 



m 

This implies that 



(C3) 
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^tYi+i)-^dt. (C4) 



By considering the fact that the exponential in the inte- 
grand is always less than or equal to 1, this formula can 
immediately be bounded above, as follows: 



jiiv) < 



2 4 



r( 



2 4> 
I 1 



rfi-i 



-B 



^(V2-l/4)-l(^_^)(i/2+l/4)-l^^ 

i± 

T{1) 



I I I i\ rC ' ^ 



2 4' 2 4 



(C5) 



where B(x,y) — T{x)r{y) /r{x + y) is the Beta func- 
tion ( p3| , Sec. 6.2). By combining this inequality with 
Eq. (pl|) we obtain the bound stated in Eq. (3.14b). 
Moreover, in the limit v <C 1, it is legitimate to replace 
the exponential factor in Eq. (C4) by 1 inside the inte- 
gral, and thus the same expression gives the asym ptotic 
form in the v -^1 regime, as quoted in Eq. (|3.13c ). 

An additional bound can be obtained for ^ > 1 by tak- 
ing into account the fact that the factor (1 — 
in the integrand is less than or equal to unity, so that 



-t.^^(//2-l/4)-l^^ = 1. (C6) 



This gives the lower bound of Eq. ( 3.14a ). When t; ^ 1, 
the same express ion pr ovides the asymptotic form for all 
values of I, Eq. ( |3.13b| ), as, in that limit, the integral is 
dominated by the region near the origin, where the factor 
(1 - t)(V2+i/4)-i is close to unity. 
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malized function tp^{k) = V"^^Sf. q concentrated at the 
zero replica sector is given by 
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